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ABSTRACT 

IRC +10216 is a circumstellar envelope around a carbon-rich evolved star which contains a large 
variety of molecules. According to interferometric observations, molecules are distributed either con¬ 
centrated around the central star or as a hollow shell with a radius of ^15". We present ALMA Cycle 0 
band 6 observations of the J = 14 — 13 rotational transition of CH3CN in IRC +10216, obtained with 
an angular resolution of 0!'76 x 0!'61. The bulk of the emission is distributed as a hollow shell located 
at just ^2" from the star, with a void of emission in the central region up to a radius of ^1". This 
spatial distribution is markedly different from those found to date in this source for other molecules. 

Our analysis indicate that methyl cyanide is not formed neither in the stellar photosphere nor far 
in the outer envelope, but at radial distances as short as 1-2", reaching a maximum abundance of 
^ 0.02 molecules cm“^ at 2" from the star. Standard chemical models of IRC +10216 predict that the 
bulk of CH3CN molecules should be present at a radius of ^ 15", where other species such as polyyne 
radicals and cyanopolyynes are observed, with an additional inner component within 1" from the star. 

The non-uniform structure of the circumstellar envelope and grain surface processes are discussed as 
possible causes of the peculiar distribution of methyl cyanide in IRC +10216. 

Subject headings: astrochemistry — line: identification — molecular processes — stars: AGB and 
post-AGB — circumstellar matter — stars: individual (IRC +10216) 


1. INTRODUCTION 

The chemical structure of the well known carbon-rich 
envelope IRC +10216 and, in general, of circumstellar en¬ 
velopes around asymptotic giant branch (AGB) stars is 
in general terms well described by an scenario in which 
stable molecules are formed in the warm and dense sur- 


round ings of the star under chemical equilibrium (Tsuji 
1973) while radicals and more exotic species are produced 
in the outer layers due to the photochemistry induced by 
the penetration of interstellar ultraviolet (UV) photons 
(|Glassgold||1996|). 

However, m recent years there have been evidences of 
a growing number of chemical aspects that do not fit into 
this general scenario. Perhaps, the most clear example is 
the detection of warm water vapour in IRC +10216 and 


other carbon-rich envelopes (iDecin et ^ 20101 Neufeld 


CJN in oxygen- 


rich envelopes ( 

Bujarrabal et al. 

1994 Bieging et al. 

2000D, NH3 in 0 

►- and C-rich enve 

opes (Ready V Ridg- 

way|l993l Menten et al.|2010), and PH3 in IRC +10216 


molecules are formed in the inner regions of the envelope 
with abundances much larger than predicted by chemical 
equilibrium. Indeed, these inner layers are complex re¬ 
gions where a variety of non-equilibrium proc esses such 
as shock s driven by the pulsation of the star (Cherchn- 
eff 2OI2I), processes related to t he condensation of dust 
grains {Gail ^ ISedlmayr||1988|), and even photochemi¬ 
cal processes driven by interstellar UV photons able to 
penet rate through the clumpy envelope (Agundez et al. 
2010), can affect the abundances of some species. Yet, 
the role of these processes on the chemistry has to be 


fully understood. 

The use of (sub-)mm interferometers able to probe the 
distribution of different molecules in the inner regions of 
circumstellar envelopes are a very promising tool to un¬ 
veil the role of the non-equilibrium processes at work 
in these inner regions. Some recent works have pre¬ 
sented interferometric observations of the inner envelope 
of IRC +10216, revealing the compact emi ssion distri¬ 
bution of some molecule s around the star (Patel et al 
20111 |Eonfrfa et al.||2014|) . The Atacama Large Millime¬ 
ter Array (ALMA) has also started to provide observa¬ 
tions of the chemical complexity in the inner envelope of 
IRC -f-10216 with a high angular reso lution (Cernicharo 


et al. 2013 
ALivA Gyc 


Velilla Prieto et al. 2015). Here we present 
eO band 6 observations of IRC +10216 with 


sub-arcsec resolution and report a peculiar emission dis¬ 
tribution for the J = 14 — 13 rotational transition of 
CH3CN. 

2. OBSERVATIONS 

The ALMA Cycle 0 observations of IRC+10216 were 
carried out in April 2012 using receiver band 6. P art of 
these observat i ons h as been previously presented by |Cer-| 
nicharo et al. (2013) and by Velilla Prieto et al. (2015). 
Here we are concerned with the J = 14 — 13 rotational 
transition of CH3CN, whose strongest K components lie 
at 257.5 GHz, within the spectral setup covering the 
255.3-260.2 GHz frequency range. 

An array of 16 antennas with baselines ranging from 
15.7 m to 384.1 m was used in three observing tracks of 
21 minutes of correlations on source. The field of view 
(primary beam) of the ALMA 12m antennas has a diam¬ 
eter of ^ 23" at 257.5 GHz. The adopted spectral reso- 





































































2 


Agiindez et al. 



Fig. 1.— ALMA uv plane coverage (left) and associated dirty beam (right) for the spectral setup covering the range 255.3-260.2 GHz. 
The primary beam of the ALMA 12m antennas at 257.5 GHz ('^ 23'') is indicated as a white circumference in the right panel. 


lution is 0.98 MHz, with a channel spacing of 0.49 MHz. 
The bright point-like source J0854+201 was observed to 
calibrate the bandpass, while the amplitude and phase 
were calibrated by observing J0854+201 and J0909+013 
every 10 and 20 minutes, respectively. The uncertainty 
in the flux calibration is estimated to be 8 %. Calibra¬ 
tion was done with the software CASAEJand further data 
processing was done using the package GILDA^ A good 
coverage of the uv plane was obtained, with 156,865 vis¬ 
ibilities (see left panel in Fig. [^. The continuum was 
subtracted by carefully selecting spectral windows free of 
line emission located in spectral regions near the CH 3 CN 
J = 14 — 13 lines. The associated dirty beam has mi¬ 
nor contributions from sidelobes (less than 10 % of the 
main beam, see right panel in Fig. [B. Image deconvolu¬ 
tion was carried out using the Ho^om clean algorithm 
with no support, i.e., no a priori structure was assumed 
for the CH 3 CN J = 14 —13 brightness distribution. This 
choice is adequate because the sidelobes of the dirty beam 
are not important and because CH 3 CN line emission in 
IRC +10216 does not have a simple structure. The re¬ 
sulting synthesized beam is 0!'76 x 0^'61 and the final rms 
per 0.49 MHz channel is 3.6 mJy beam“^. 

We carried out observations with the IRAM 30m tele¬ 
scope in June 2014 to estimate the degree of flux that 
could have been filtered out by ALMA. We used the 
EMIR E230 receiver connected to a fast Eourier Trans¬ 
form spectrometer providing a spectral resolution of 0.2 
MHz (data were later on resampled to the channel spac¬ 
ing of 0.49 MHz of the ALMA data). We performed a 
scan over a region 44" x 44" during 18 h of on source in¬ 
tegration time using the ‘onthefly’ observing mode. The 
half power beam width (HPBW) of the IRAM 30m tele¬ 
scope at 257.5 GHz is 9"4. System temperatures were in 
the range 300-350 K, resulting in a rms of ^0.02 K 
(0.17 Jy) per 0.49 MHz channel in the central region 
of the map. The spectrum obtained with the IRAM 
30m telescope at the position of the star is compared 

^ See http://casa.nrao.edu 

^ See http://www.iram.fr/IRAMFR/GILDAS 


in Eig. with the ALMA spectrum integrated over a re¬ 
gion 9M X 9"4, i.e., similar to the HPBW of the 30m 
telescope. It is seen that the flux filtered out by ALMA 
is negligible at the line edges of the different K lines, 
which correspond to the terminal expansion velocity of 
the envelope and therefore show a rather compact distri¬ 
bution, although it can be as high as 30-40 % in some 
channels located around the center of the K lines, where 
the emission is more extended. 


3. RESULTS 

Some conclusions can be drawn after inspecting the 
brightness distribution in the different channels cover¬ 
ing the most intense K components {K = 0 — 3) of the 
GH 3 GN J = 14 — 13 transition. Hereafter, we refer to 
^LSR ~ ysysi he., the local standard rest velocity cor¬ 
rected for the systemic velocity of the source, as simply 
V. In the case of IRG +10216, the systemic velocity Vgys 
is —26.5 km s“^ and the ter minal expansion veloci ty of 
the envelope is 14.5 km s“^ (Gernicharo et al.|| 2000 ). At 


velocities close to the terminal expansion velocity, the 



Frequency (GHz) 

Fig. 2.— Spectrum of IRG+10216 in the region of the GH3GN 
J = 14 — 13 transition as observed with the IRAM 30m telescope 
and with ALMA (integrated over a field of view of similar size to 
the HPBW of the 30m telescope). 
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Fig. 3.— Observed brightness distributions for the AT = 0,1, 3 components of CH 3 CN J = lA — 13, averaged over a velocity range of width 
2.8 km s“^ centered at Vlsr = Vsys- Since each line has a width of ~29 km s“^, the K=0 and K=1 components overlap and therefore the 
Jx = 14o — 13o map contains a contribution from the transition K = 1 line at a velocity V = —5.8 km s“^ while the Jk = 14i — 13i map 
is also contaminated by emission from the AT = 0 component at y = +5.8 km s“^. The Jk = 143 — 133 map (also shown after merging 
ALMA and 30m short-spacing data; bottom right panel) is free of contamination and therefore traces better the spatial distribution of 
CH 3 CN in the plane of the sky. The rms per 2.8 km s“^ channel is 1.6 mjy beam“^ and contour levels are shown at 5, 10, 15, and 20 times 
the rms. The maps are centered on the star, indicated by a cross and with coordinates J2000.0 R.A.=09^47"^57?445, Dec.=13°16h43''86. 
The size and shape of the synthesized beam (0''76 X 0''61) is shown in the bottom left corner of each panel. 


emission appears compact and centered on the star while 
at velocities around Vlsr = Vsys the emission is dis¬ 
tributed as a ring around the central position. This is the 
typical pattern of a molecule which is not concentrated 
around the central star but distributed in a hollow shell 
at a certain distance from the star. We must however 
note that in the case of CH 3 CN J = 14 — 13, this pat¬ 
tern is less clearly appreciated for the K = {) and K = 1 
components because these two lines overlap (their rest 
frequencies are closer than the line width of 29 km s“^), 
and therefore some of the channels have contributions 
from both K components at different velocities. 

In Fig. we show the brightness distribution of the 
strongest K components {K = 0,1,3) of the J = 14 — 13 
transition of CH 3 CN at Vlsr = Vsys- We first focus on 
the maps based on ALMA data without including short¬ 
spacing (top and bottom left panels in Fig. [^. It is seen 
that the maximum of the emission appears as a ring with 
a radius of ^ 2" and that there is a hole within ^ 1" from 


the position of the star. The map corresponding to the 
Jk = 14o — 13o transition contains also emission from the 
K = 1 component at R = —5.8 km s“^, which makes the 
ring to broaden inwards. The same occurs in the map 
of the Jk = 14i — 13i transition, which is contaminated 
by the K = ^ component at R = +5.8 km s“^. In the 
case of the K = 2 and K = 3 lines, emission at Vlsr = 
Vsys is free of blending with other K components. In 
particular, the map of the Jk = 143 “ 133 transition, 
which is more intense than the Jk = 142 ~ 132 transition, 
allows to see clearly the ring structure. It is also worth 
to note that a second ring of emission located at ^ 4" 
from the star is barely apparent at a few a in the Jk = 
14i — 13i map, suggesting a spiral structure which has 
been also observed with ALMA in other molecular lines 
(Cernicharo et al. in prepa ration) and at large sca les in 
the J = 2 — 1 line of CO (| Cernicharo et al. 2015|), and 
that has been interpreted in terms of a binary system. 

To evaluate whether the brightness distributions of the 
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Fig. 4.— Brightness distribution of CH 3 CN Jx=14o-13o, aver¬ 
aged over a velocity range V from —2.5 to +2.5 km s“^, as would be 
seen with a synthesized beam of 5^9 x 3^0 (shown in the bottom 
left corner). To be compared with the brightn ess distribution of 
the CH 3 CN Jic=16o-15o transition observed by|Patel et al.|(|2011|) 
with SMA (line 9 in their Figure 5). Note that with tins angu¬ 
lar resolution the inner hole of emission seen with ALMA is not 
resolved. 


different K lines at Vlsr = Vsys could be significantly 
affected by the ALMA flux loss, we have combined the 
ALMA data with the IRAM 30m short-spacing data to 
include visibilities at shorter baselines than the shortest 
one in the ALMA data (15.7 m). We however note that, 
given the modest sensitivity of the IRAM 30m data and 
the fact that the CH 3 CN emission is restricted to a region 
much smaller than the primary beam of ALMA 23"), 
we do not expect great changes in the emission struc¬ 
ture. As an example, we compare in the bottom panels 
of Fig.j^the brightness distribution of the Jk = 143 — 133 
line at Vlsr = Vsys^ as obtained using ALMA data alone 
and after combining ALMA and 30m short-spacing data. 
It is seen that, appart from a slight recovery of flux at ex¬ 
tended scales, the emission structure remains essentially 
unaltered. 

Higher J lines of CH 3 CN were observed in IRC +10216 


with t he Submillimeter Array (SMA) by Patel et al. 
( 2011 ) in their spectral survey covering the frequency 
range 293.9-354.8 GHz. These authors find a compact 
brightness distribution for the CH 3 CN lines, observed at 
angular resolutions ranging from 2 " to 6 ", i.e., noticeably 
poorer than the angular resolution of the ALMA obser¬ 
vations presented here. Moreover, these authors assign 
various emission lines to rotational transitions of CH 3 CN 
in the z/g = 1 vibrational state on the basis of a good 
agreement between the observed and laboratory rest fre¬ 
quencies. These results may lead to the conclusion that 
in IRC +10216, CH 3 CN is present in the warm surround- 



that should be more intense than those identified are 
missing. For example, in the case of the J = 16 — 15 


TABLE 1 

Vibrational modes of CH 3 CN 


Vibrational mode 

Sym.“ 

Freq.^ 

(cm“^) 

A" 

(s-i) 

ui CH 3 symmetric stretching 

A 

2954 

3.1 

1^2 CN stretching 

A 

2268 

1.1 

1^3 CH 3 symmetric deformation 

A 

1385.2 

- 

1^4 CC stretching 

A 

920.290 

0.21 

^'5 CH 3 antisymmetric stretching 

E 

3009 

I.O 

uq CH 3 antisymmetric deformation 

E 

1449.7 

5.3 

CH 3 rocking 

E 

I04I.855 

0.50 

Us CCN bending 

E 

365.024 

0.016 


Note. — Table built from the compilation by J. Crovisier at 
http://www.lesia. obspm.fr/perso/iac ques-crovisier/basemole 
and from Table 1 of|Miiller et al.| (|2015|). 

“ Symmetry of first excited state ot eacJi mode {vi = 1). 

^ Frequency and Einstein coefficient for spontaneous emission of 
the fundamental band of each mode (i/j = 1 —)■ 0 ). 


transition in the ug = 1 state, only the K = 11,12,13 
components are identified, but not the K = 0,1,2,3 
lines, which should be far more intense. Therefore, the 
lines assigned to CH 3 CN z/g = 1 must arise from other 
carriers. On the other hand, as concerns the rotational 
lines of CH 3 CN in the ground vibrational state observed 
with SMA, the compact emission distribution observed 
is probably a consequence of the limited angular reso¬ 
lution, which does not allow to resolve the inner hole 
seen by ALMA. This is illustrated in Fig. where we 
have applied a tapering procedure to the ALMA data to 
suppress long baselines and simulate what would be seen 
with a synthesized beam of 5"9 x 3"0, as i n the SMA 


obser vations of, e.g., the Jk = I 60 — 15o (see Patel et al. 


2011). It is seen that with such limited angular resolution 
the inner hole in the emission is not resolved. Moreover, 
since the CH 3 CN lines observed by SMA are higher in J 
than the one observed with ALMA, their emission could 
be shifted inwards due to excitation requirements, and 
thus a higher angular resolution would be needed to re¬ 
solve the inner hole. We thus conclude that the hollow 
shell structure seen by ALMA in CH 3 CN J = 14 — 13 
emission is not in conflict with the compact distribution 
observed for higher J transitions with SMA. 


4. DISCUSSION 

According to the interferometric maps carried out 
during the past 20 years, the molecular emission in 
IRC +10216 appears either concentrated around the cen¬ 
tral star, as in the case of HCN, CS, SiO, SiS, and NaCN, 
or distributed in a hollow shell located at a radius of 


^ 15" from the star, as occurs for HC 5 N, H NC, MgNC, 


and the radicals CN, C 2 H, C 3 H, and C 4 H (|Lucas et al. 

19951 

Dayal & Bieging||1995l IGuelin et aL||1997[ |Fbntria 

et al. 

20i4| Velilla Prieto et al. 2015). However, in the 

r t M i /mi xi__II _ X- 1 _ X - ] ^x- lir// _ 


case of GHgCr^, the shell is not located at ^ 15" from 
the star but at just ^ 2 ". 


4.1. Excitation and radiative transfer model 

In order to use the ALMA observations to put con¬ 
straints on the abundance, distribution, and excitation 
of CH 3 CN in the envelope of IRC +10216 we have car¬ 
ried out excitation and radiative transfer calculations 
based on a multi-shell LV G formalism. We adopt ed the 
model of the envelope of Agiindez et al. ( 2012 ), with 
the downward revision of the density of particles de- 
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Fig. 5.— Radial abundance profiles of CH 3 CN in IRC+10216 
are given as fractional abundance relative to H 2 in the upper panel 
and as volume density, scaled in some cases, in the lower panel. 
Solid lines represent the abundance profile that best reproduces 
the observations, dashed lines the result of a standard chemical 
model of the envelope, and dotted lines refer to a photochemical 
model of a clumpy envelope. The radius is given in arcseconds for 
an assumed distance of 130 pc. 


rived by Cernicharo et al. (2013|) for the innermost layers 
(< 5R*). T'he adopted mass loss rate and dis tance are 
2 X 10“^ Mq yr“^ and 130 pc (see Table 9 of Agiindez 


|et al. 1 ( 2^2 for other parameters). As concerns the spec- 
troscopy of CH 3 CN we consider separately the A and 
E species with the A:E statistical ratio of 1:1 and com¬ 
pute the en ergy levels from the spe ctro scopic constants 


reported by 


zarini 


(2006 
moment of 


Simeckova et al. 


(2004) and 


Cazzoli & Puz- 


2008) and the line strengths from the dipole 


.92197 ± 0.00013 D measured by [Gadhi et 


al. (1995). Infrared pumping to excited vibrational states 


may also play a role in the excitation of CH 3 CN. This 
molecule has eight vibrational modes, whose fundamen¬ 
tal bands lie at wavelengths from 3.3 fim to 27 /im (see 
Table [^, and the radiation field in the highly reddened 
envelope of IRC +10216 is intense in the mid-infr ared, 
especially around 10 /im (Cernicharo et al. 1999). A 
detailed analysis of the complete pumping scheme of 
CH 3 CN is beyond the scope of this article. However, 
to have an idea of the effect of infrared pumping we in¬ 
cluded the first excited state of the vibrational mode z/g, 
which can be pumped by 6.9 /im photons. Other vibra¬ 
tional bands could also be important, e.g., the funda¬ 
mental band of z/ 3 , whose strength is unfortunately not 
available, and overtones and hot bands of low-frequency 


Fig. 6 .— Radial brightness distribution of the K = 0 to K = 3 
components of CH 3 CN J = 14 — 13, averaged over a velocity range 
of width 2.8 km s“^ centered at Vlsr = Vsys- Note that the 
K = 0 and K = 1 components contain also a contribution from 
K = 1 ai V = —5.8 km s“^ and from K = 0 at V = +5.8 km s“^, 
respectively. Solid lines correspond to the ALMA data combined 
with the IRAM 30m short-spacing data and dashed lines to the 
results of the model that best reproduces the observations (black 
and red dashed lines are almost indistinguishable). 


modes such as ug. Vibrational modes such as z^i, 1 ^ 2 , or 
z ^5 have also strong fundamental bands, although they 
lie at shorter wavelengths where the radiation field in 
IRC +10216 is less intense. The spectroscop ic constants 
of the z /6 = 1 state have been taken from [Paso et al.| 
(1994). The strength o f the vibrational band =0^1 
has been measured by Cerceau et al. (1985). We con¬ 
sidered rotational levels within the ground and 1 /^ = 1 
vibrational states up to J = 30. We adopted the rate co¬ 
efficients of pure rotational excitation thr ough inelasti c 
collisions with H 2 from the calculations by Green (1986), 
which extend up to 140 K, and extrapolated at higher 
temperatures assuming that the rate coefficients scale 
with the square root of temperature. Collisional exci¬ 
tation for ro-vibrational transitions was neglected. 

As radial abundance profile for CH 3 CN we used an 
expression of Gaussian type which takes the form 


where / is the fractional abundance with respect to H 2 , 
ro is the radius where / peaks, and Alogr is a measure 
of the radial extent. The radial abundance profile which 
best reproduces the J = 14 — 13 transition observed with 
ALMA is given by Eq. ^ with ro = 1.2 x 10^^ cm, 
/(fo) = 5 X 10“^, and Alogr = 0.45 for r < ro and 
0.3 for r > ro (see solid lines in Eig. It is seen that 
the fractional abundance relative to IC reaches a max¬ 
imum value of 5 X 10“^ at 6" (see solid line in upper 
panel of Eig. [^. However, because of the increase in the 
gas density with decreasing radius, if the abundance is 
expressed as number of CH 3 CN molecules per unit vol¬ 
ume, the peak value is reached well before, at about 2" 
(see solid line in lower panel of Eig. [^. 

The radial brightness distributions calculated for the 
various K components of the J = 14 — 13 transition are 
confronted with the observed ones in Eig. According 
to our calculations, the excitation of the rotational levels 
involved in the J = 14 — 13 transition is dominated by in¬ 
elastic collisions with H 2 , although infrared pumping to 
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the = 1 state does also play a role as it induces a slight 
increase in the line intensity of about 10 % and shifts the 
emission to slightly outer radii. Including other vibra- 
tionally excited states could enhance the importance of 
infrared pumping. We find that as the quantum number 
K increases, and so does the upper level energy of the 
corresponding transition, the excitation requirements be¬ 
come more demanding and the maximum of brightness 
shifts to shorter radii, where the higher gas density and 
infrared radiation field favor a more efficient excitation. 
This effect is however noticeable for high—K transitions. 
In the case of the K = 0 to K = 3 components, all them 
show a brightness peak at about 2 ", the same radius at 
which the volume density of CH 3 CN molecules reaches 
its maximum. It is also worth noting that the observed 
radial brightness distribution of three out of the four K 
components shown in Fig. show a more or less marked 
double-humped shape, which probably arises from the 
non-uniform character of the CH 3 CN shell, with sub¬ 
shells at smaller scales, spiral-like structures, etc. Our 
model of a smooth hollow spherical shell does not aim 
at reproducing these fine details, but the global shape of 
the observed brightness distribution. 

In any case, it is clear that the CH 3 CN shell is located 
in IRC +10216 at a much shorter radius than the molec¬ 
ular shell at ^15", where the radicals CN, C 2 H, C 3 H, 
and C 4 H, and molecules su ch as HNC and HC 5 N are 
observed ( Gnelin et al.||1997 ). 


4.2. Chemical model 

Standard chemical models of IRC +10216 cannot re¬ 
produce the ALMA observations in th at they predict 
that CH 3 CN forms at too large radii (Agiindez et al. 
20081 |Li et al.||2014|) . Here we have constructed a stan¬ 
dard chemical model using a gas phase reactio n network 
built up from the databases UMIS T and KID A (McElroy 
et al.||2013 |Wakelam et al. 2015) and the literature on 
gas phase memical kinetics. We have chosen an initial 
radius of 2 x 10 ^^ cm, which corresponds to O^'l for an 
assumed distance of 130 pc, to investigate the possible 
formation of CH 3 CN in the inner regions of the envelope. 
The initial abundance s injected into the exp anding en¬ 
velope are taken from Agiindez et al. ( 2012 ) and N 2 is 
assumed to have an abundance of 4 x 10“^ relative to 
H2. 

According to the chemical model, the formation of 
methyl c yanide throughout th e envelope of IRC+10216 
(see also Agiindez et al.||2008|) involves the formation of 
the precursor ion CH 3 CJVH+ through the radiative asso¬ 
ciation reaction 


CH+ + HCN ^ CH 3 CNH+ + hv, 


( 2 ) 


followed by the dissociative recombination of CH 3 CNH+ 
with electrons. Among the reagents in reaction ( 0 , HCN 
is fairly abundant from t he innermost layers out to the 
photodissoc iation region ( Fonfrfa et al.|2QQ8 Cernicharo 
et al. 2011 ) while the ion CHj appears from the first 
stages of the expansion driven by cosmic rays although 
the bulk of formation is induced by photochemistry in the 
outer layers. The calculated abundance of CH 3 CN rela¬ 
tive to H 2 is shown as a function of radius in the upper 
panel of Fig. (dashed line). Note the early appearance 
driven by cosmic rays at radii shorter than ^ 10 ^^ cm and 


the formation bump induced farther out by photochem¬ 
istry. As previous chemical models, ours also puts most 
of CH 3 CN in the outer envelope, in a shell located at 
about 15" from the star. It is however interesting to note 
that the non-negligible fractional abundance of CH 3 CN 
driven by cosmic rays from the inner layers together with 
the high gas densities prevailing in these regions lead to a 
distribution of CH 3 CN molecules which, appart from the 
component located at 15", has also an inner component 
within 1" from the star. This is more clearly appreciated 
when expressing the abundance in absolute rather than 
relative terms (dashed line in lower panel of Fig.[^. The 
finding of the inner component of methyl cyanide in our 
chemical model is interesting because it appears closer 
to the region where this molecule is observed than the 
outer shell at 15". However, this inner component would 
appear compact rather than showing an inner hole, for 
angular resolutions poorer than ^ 0 " 1 , as is the case of 
our ALMA data. 

Our chemical model underestimates the abundance of 
CH 3 CN by more than two orders of magnitude and pre¬ 
dict a radial distribution markedly different from that 
observed. There are however uncertainties in the chem¬ 
ical model which are worth to discuss. Concerning re¬ 
action rate constants, that of the dissociative recombi¬ 
nation of CH 3 CNH+ with electrons is known from mea¬ 


sur ement s involving the fully deuterated ion (Vigren et 


al. 2008), although the branching ratios of the difier- 


ent channels, inclu ding that yielding C H 3 CN, are based 
on simple guesses_([Vigren et al. 2009). The rate con¬ 
stant of reaction (0 is not well constrained. Low pres¬ 
sure ex periments sug gest a value of 2 xl 0 “^^ cm^ s“^ at 
300 K ( Anicichj|1993) while statistical cal culations point 
to a value 45 times nigher (Herbst||1985). In the chem¬ 
ical model we adopt the lower experimental value. If 
the higher theoretical value is adopted, the calculated 
abundance of CH 3 CN experiences a ten-fold enhance¬ 
ment, although the shape of the radial profile remains 
similar. Another important parameter that affects the 
abundance of methyl cyanide, especially in the regions 
inner to 10 ^^ cm, is the cosmic-ray ionization rate of 
H 2 , which in the chemical model is set to a value of 
1.2 X 10“^^ s“^. If we adopt a higher rate, the abun¬ 
dance of CH 3 CN increases at short radii. That is, the 
inner component gains importance with respect to the 
outer one. In any case, the ionization rate due to cosmic 
rays in IRC +10216 is relatively well constrain ed by ob 


servations of HCO^ at millimeter wavelen gths (Agiindez 


& Cernicharo 2006 Pulliam et al. 2011), and therefore 


^ ^ s ^ 


cannot be much higher than a few 10 
As long as a standard chemical model fails to explain 
the distribution of CH 3 CN in IRC +10216, it seems clear 
that such model misses either some key chemical process 
or an adequate description of the physical structure of the 
envelope at arcsecond scales. There are various aspects 
which are worth to discuss. It is known that the enve¬ 
lope IRC +10216 is not completely smooth, but t hat it 


consists of shells and arcs which cross each ot her (Mau- 

le ef- 


ron & Huggins 
feet of these she 


2000 Cernicharo et al. 2015) 


Is on the radial distribution o: 


T 

■ molecules 


IRC +10216 has been studied by [Cordiner fc Mill^ 


(2009), although the case of CH 3 CN was not specihcally 
discussed by these authors. A non-uniform radial density 
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profile could modulate the early appearance of CH 3 CN 
in the envelope, making the inner component to appear 
as a small hollow shell rather than a compact structure. 
To investigate this possibility it is necessary to have a 
precise knowledge of how the gas is distributed in the 
inner circumstellar layers. Another aspect of having a 
non-uniform distribution of matter in the circumstellar 
envelope is the clumpy nature, which may allow interstel¬ 
lar UV photons to penetrate in the inner layers, allowing 
to form some molecules as a result of p hotochemistry 


(Decin et al. 2010 Agundez et al. 


2010). In Fig. ^ we 


show the expected distribution of CH 3 CJN under such sce¬ 
nario a dopting the same p arametric treatment of dumpi¬ 
ness of Decin et al. (2010). It is seen that the abundance 
of CH 3 CW indeed increases in the inner layers (dotted 
lines in Fig. [^, although the distribution of molecules 
would appear compact around the star rather than as 
a hollow shell. It remains to be seen whether a dif¬ 
ferent scenario of clumpy envelope or a combination of 
dumpiness and shells could result in an abundance dis¬ 
tribution of CH 3 CN compatible with that observed by 
ALMA. Alternative formation mechanisms such as grain 
surface processes could also play a role in the synthesis 
of CH 3 CN. Such possibility, although difficult to prove, 
is an interesting one because CH 3 CN appears soon after 
the dust condensation region in IRC +10216. Additional 
dues to understand the puzzling distribution of CH 3 CN 
in IRC +10216 could come from further ALMA observa¬ 
tions of other species showing similar distributions. In 
particular, it would be interesting to see whether the 
chemically related species CH 2 CN shows a similar spa¬ 
tial distribution, although the lines of this radical are 


perhaps too faint to be mapped ( Agundez et al.|[2008 ). 


5. CONCLUSIONS 

We used ALMA to observe the J = 14 — 13 rota¬ 
tional transition of CH 3 CN, lying at 257.5 GHz, with 
an angular resolution of 0!'76 x 0^61 in the carbon star 
envelope IRC +10216. The bulk of the emission is dis¬ 
tributed as a hollow shell ^ 1" wide located at ^ 2" 
from the star. The observations, together with an excita¬ 
tion and radiative transfer analysis, indicate that methyl 
cyanide is not formed neither in the stellar photosphere 
nor far in the outer envelope, but at radial distances 
as short as 1 — 2 ", reaching a maximum abundance of 
^ 0.02 molecules cm“^ at 2" from the star. This fact is 
in conflict with standard chemical models of IRC +10216, 
which predict that the bulk of CH 3 CN molecules should 
be present at a radius of ^ 15" with an additional inner 
component within 1 " from the star. It is not yet clear 
which mechanism is at the origin of the peculiar spatial 
distribution of methyl cyanide in IRC +10216. The non- 
uniform structure of the envelope, with shells, arcs, and 
clumps, and grain surface processes may play an impor¬ 
tant role on the distribution of CH 3 CN. Further ALMA 
observations of other molecules showing spatial distribu¬ 
tions similar to that of CH 3 CN should bring light on the 
puzzling distribution of methyl cyanide in IRC +10216. 
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